EMPG-04-02 



Auxiliary matrices for the six-vertex model 
and the algebraic Bethe ansatz 

Christian KorfT 

School of Mathematics, University of Edinburgh 
Mayfield Road, Edinburgh EH9 3JZ, UK 



Abstract 

We connect two alternative concepts of solving integrable models, Baxter's 
method of auxiliary matrices (or Q-operators) and the algebraic Bethe ansatz. 
The main steps of the calculation are performed in a general setting and a formula 
for the Bethe eigenvalues of the Q-operator is derived. A proof is given for states 
which contain up to three Bethe roots. Further evidence is provided by relating 
the findings to the six- vertex fusion hierarchy. For the XXZ spin-chain we analyze 
the cases when the deformation parameter of the underlying quantum group is 
evaluated both at and away from a root of unity. 
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1 Introduction 



This paper is a continuation of two previous works Ej on the six-vertex model and 
the associated XXZ Heisenberg spin-chain at roots of unity. That is, we consider the 
integrable model defined via the Hamiltonian 

M + x 

H = a m a m+l + Vm+1 H 4 ^m^m+l; ^M+l = A ±2 of , CT^ +1 = of . (1) 

m=l 

Here {a x , a y , a z } are the Pauli matrices with <j ± = (cr 1 ± ia y )/2 and A G C fixes the 
boundary conditions. The anisotropy parameter in front of the third term in (JTJ is fixed 
in terms of the complex variable q which in general will be of modulus one. Of particular 
interest to our discussion is the case when q is a primitive root of unity, i.e. q — 1 for 
some integer N > 2. ( We exclude here the cases N = 1,2 which are related to the XXX 
model.) At these particular values the above Hamiltonian as well as the associated six- 
vertex transfer matrix exhibit extra degeneracies in their spectra which are linked to an 
underlying infinite-dimensional non-abelian symmetry algebra. At periodic boundary 
conditions, A = 1, and in the commensurate sectors 2S Z = OmodA this symmetry 
algebra is isomorphic to the loop algebra of sl2 ■ This has been first established by 
Deguchi, Fabricius and McCoy in [3], where additional results for N = 3,4 outside the 
commensurate sectors can be found. When A ^ 1 the symmetry algebra will in general 
reduce to the upper or lower Borel subalgebra [HE]- A special case is obtained when 
the boundary conditions are tuned to A = q ±s * . Then the symmetry algebra and the 
explicit form of its generators are known for all spin sectors and integers N jS]. 

The presence of an non-abelian symmetry algebra is of interest as it now allows one 
to connect the Bethe ansatz and representation theory. Note that while this idea is 
not new, the mentioned symmetries are infinite-dimensional and exist at finite length 
of the chain, M < 00. The combination of these two distinct features distinguishes the 
present discussion from previously considered cases in the literature at open boundary 
conditions jH] or at infinite volume jjj. 

An important first step in connecting the Bethe ansatz with the representation theory 
of the aforementioned symmetry algebras is to find an efficient way to analyze the 
structure of the degenerate eigenspaces of the Hamiltonian as well as the transfer matrix. 
This can be achieved by using Baxter's concept of auxiliary matrices [S], also known as Q- 
operators, which satisfy certain operator functional equations with the transfer matrix. 
The concept of auxiliary matrices has primarily received attention in the context with 
the eight- vertex model [HJ HH1 EEU U2] ; see [T3J Ej for a recent discussion corresponding 
to the root of unity case. 

In connection with the six-vertex model the subject has obtained fresh impetus from 
new methods of constructing such Q-operators which provide an alternative to Baxter's 
procedure described in e.g. [TS]. These new methods are based on the representation 
theory of quantum groups [THl HZ1 UH EE EE] ■ The latter method will be of importance to 
us as the auxiliary matrices constructed in this manner, see PUI2|j have been shown for 
several examples to yield information on the irreducible representations of the symmetry 
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algebra at roots of unity. Several facts about the spectrum of the auxiliary matrices in 
[2] have so far only been rigorously proven for N = 3 and conjectured to hold true for 
N > 3 employing numerical calculations. 

The purpose of this article is to lend further support to the earlier conjectures and to 
extend the discussion from auxiliary matrices with periodic boundary conditions (A = 1) 
to quasi-periodic ones (A ^ 1) in order to accommodate the findings in Sj. To this end 
we take a broader point of view and consider the connection between the algebraic Bethe 
ansatz [20] and auxiliary matrices in a general setting and at generic values of q. 

Away from a root of unity this will enable us to calculate the spectrum of the 
auxiliary matrices constructed in (TH] and resolve certain convergence problems due to 
an infinite-dimensional auxiliary space. We will also make contact with the discussion 
in [TU. 

At roots of unity |2| the Bethe ansatz analysis will not yield the complete set of 
eigenvalues for the Q-operators, as at most the highest weight state in each degenerate 
eigenspace of the transfer matrix ought to be a proper Bethe state, i.e. a state parame- 
terized by finite solutions of the Bethe ansatz equations. Detailed explanations will be 
given in the text. In addition to the comparison with the algebraic Bethe ansatz we will 
also make contact with the fusion hierarchy of the six-vertex model. The latter provides 
an infinite series of higher-spin transfer matrices which can be successively generated 
through a functional equation. At roots of unity this series truncates and we will show 
how the auxiliary matrices are related to the fusion matrices, similar to the discussion in 
[Tl] for the eight-vertex model. This will provide additional evidence for the spectrum 
of the auxiliary matrices constructed in QjJ |2] . 

The paper is organized as follows. In Section 2 we introduce the monodromy matrices 
of the six- vertex transfer matrix and the Q-operators constructed in [I9J and 
Afterwards we derive the commutation relations between the auxiliary matrices and the 
Yang-Baxter algebra. 

This sets the stage to compute the action of the auxiliary matrices on Bethe states in 
Section 3. We will find that they are eigenstates of the auxiliary matrices and compute 
the corresponding eigenvalues. This is done without specifying the quantum space. The 
computations are quite lengthy and cumbersome whence we will only give a proof for 
Bethe states which contain up to three Bethe roots. For the general case of an arbitrary 
number of Bethe roots we formulate a conjecture. 

In Section 4 we specialize our findings to the XXZ spin-chain. While we discuss both 
q being a root of unity and q not being a root of unity, the former case is considered 
in more detail in light of the aforementioned symmetries. We check the conjectured 
formula for the eigenvalues of the Q-operators for consistency by inserting them into 
the respective functional equation with the transfer matrix. 

In Section 5 we present further support for the conjecture regarding the spectrum 
of the auxiliary matrices by showing that the eigenvalues are also consistent with the 
fusion hierarchy of the six-vertex model. The fusion hierarchy is solved explicitly in 
terms of Bethe roots and we derive its "truncation" at roots of unity. 

Section 6 contains the conclusions. 
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2 Monodromy matrices 



The first step in the application of the algebraic Bethe ansatz is the definition of the 
monodromy matrices. The monodromy matrices act on a tensor product of two spaces, 
'Hq <8> Ti.] the first is called the auxiliary apce and the second the quantum space. De- 
pending on the choice of the auxiliary space we obtain the monodromy matrices for 
the Q-operator and the six-vertex transfer matrix. Since the monodromy matrices are 
subject to the Yang-Baxter equation they can be constructed in the framework of quan- 
tum groups. Let {e^, ft, q hi } i=01 be the Chevalley-Serre generators of the quantum loop 
algebra U q (sl 2 ) subject to the relations 

q hi ejq~ hl = q Aij ej, q hl fjq~ hl = q~ Alj fj, q hl q hj = q hj q h \ i,j = 0, 1, (2) 

where the Cartan matrix reads 



A 



2 -2 
-2 2 



We will only be dealing with representations where the central charge of the affine 
extension is zero, whence we set 

h = hi — —h . (3) 
In addition, for i ^ j the Chevalley-Serre relations hold, 

e i e i - Wi^i + \^\ q ^A - = o, 

fifi-® q f?fifi + [Wifi-fifi =0. (4) 

The quantum algebra can be made into a Hopf algebra, its most important property 
being that of a coproduct which we choose to be (i = 0, 1) 

A(e i ) = l®e i + q tH ®e i , A(/j) = ft <g> q~ hi + 1 ® f h A(q hi ) = q hi <g> q hi . (5) 

The opposite coproduct A op is obtained by permuting the two factors. If the deformation 
parameter q is considered to be an abstract indeterminate there exists the universal R- 
matrix intertwining these two coproduct structures 

RA(i) = A op (x) R, x g U q (sl 2 ), R e U q (b+) ® U q (b.) . (6) 

Here U q {b±) denote the upper and lower Borel subalgebra, respectively. We now define 
the six-vertex monodromy matrix by setting 

T(z) = (tt^ ® 7i n )\ hm R e End C 2 ® H, (7) 

with tt^P : U q {sl2) — > EndC 2 being the fundamental evaluation representation given in 
terms of Pauli matrices, 

Tr^eo) = za~, *P(fo) = z- 1 a\ n^(q h °) = q~° z , 

7rf)( ei ) = a\ ^\h) = a-, ^\q^) = f . (8) 
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We will denote the basis in the corresponding representation space C 2 by {|0) , |1)}. 
The representation tt-h in the second factor determines the quantum space of our theory. 
At the moment we leave it unspecified in order to emphasize the general nature of the 
following discussion. Later we will set ft = (C 2 )® M and 7r^ = <S^=i 7r [ m order to 
obtain the inhomogeneous XXZ spin-chain. The factor in front of the universal i?-matrix 
involving A fixes the boundary conditions. 

In the context of the algebraic Bethe ansatz it is customary to decompose the mon- 
odromy matrix w.r.t. the auxiliary space C 2 , i.e. one introduces the following elements 
of End ft, 

A=(0|T|0) c2 , 5=(0|T|1) C2 , C = (1|T|0) C2 , £> = (1|T|1> (? . (9) 

Here the subscript C 2 indicates that the matrix elements are taken w.r.t. the first 
factor in (JZj). These elements obey the familiar commutation relations of the six- vertex 
Yang-Baxter algebra which are deduced from the relation 

R 12 (w/z)T 1 (w)T 2 (z) = T 2 (z)T 1 (w)R 12 (w/z), (10) 

with 

R(z, q) = 2±*1 ®l + ^a z ®a z + ca + ®a~ + c'a~ <g> a + (11) 

denoting the six-vertex i?-matrix. The parameterization of the Boltzmann weights is 
chosen in accordance with © and (JEJ), 

q, c = - -, c = cz . (12) 



1 — zq z 1 — zq A 

We next define the monodromy matrix from which we will obtain below the auxiliary 
matrix or Q-operator. We follow an analogous procedure to that just given with the 
difference that we now change to a higher dimensional auxiliary space fto belonging to a 
suitably chosen representation n' w : U q {b + ) — > Endft - The auxiliary space TLq will turn 
out to be finite-dimensional at q N = 1 and infinite-dimensional when q is not a root of 
unity. The monodromy matrix is now 

Q(w) = ( n ' w <g> n n )X mi R e End H <8> H . (13) 

At roots of unity there arises a technical subtlety. The universal R-matrix need 
not exist as whether one can construct an intertwiner depends on the precise nature of 
the representations tt' w ,it-h. This is due to the enlarged centre of the quantum group 
at roots of unity, see e.g. 21. However, we will assume that the representations are 
always chosen such that this is the case. Then the definition ()13|) has to be replaced by 
this intertwiner which is to be explicitly constructed. (In the case of the monodromy 
matrix ((7j) one can first evaluate it away from a root of unity and then safely take the 
root of unity limit.) 

Similar to the six-vertex monodromy matrix we may also decompose (JT5j) over the 
auxiliary space, 

Q = (Qy) with Qij := (i\Q,\j) Ho e End ft, (14) 
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where the states are labelled by integers lying in a finite interval, < i,j < N' — 1, at 
roots of unity q N = ±1 or in an infinite interval i,j £ Z when q N ^ ±1. Here and in 
the following we set N' = N if the order of the root of unity is odd and N' = N/2 when 
it is even. 

We will specify the respective auxiliary spaces momentarily. Before so doing we 
introduce a third operator L e End7io ® C 2 which intertwines the tensor product 
tt' w <8> 7r£ 1 ' ) and satisfies together with the monodromy matrices (J7|) and (JTHj) the Yang- 
Baxter equation, i.e. 

L 12 (w/z)Q 1 (w)T 2 (z) = T 2 (z)Q 1 (w)L 12 (w/z) . (15) 
Here we have assumed that we are only allowing for representations ir' w in (II 3(1 for which 

[L(w),ir' w (h) <g> 1 + = . (16) 

This restriction is necessary in order to accommodate the quasi-periodic boundary con- 
ditions when A / 1. Again, we decompose the L-operator, this time over the second 
factor, the two-dimensional space associated with the evaluation representation ir^, 

L=(" gj, a,p,~f,6 eEndHo ■ (17) 

We shall now explicitly specify the evaluation representations ir f w at and away from a root 
of unity as well as the matrix elements of the L-operator within these representations. 



2.1 The auxiliary space when q = 1 

We adopt the conventions in j2] and define the following evaluation representation n' w = 
7T^, of the full affine quantum algebra U q (sl 2 ) at q N — 1 or equivalently q N = ±1. Let 
< n < N' - 1, n e C and set 

<(fi)\n) = \n+l), <(eo) =«;<(/!), 

<(ei)|n> = ^ + ^ l ^ q ^~y' lq ' 2n \n~l), <(/o) = uT 1 

|n> = //-V^ln), <(9 h0 )=<(g- fcl ). (18) 
The non-vanishing matrix elements of the L-operator then read P 

a n = (n\a\n) = (w/z) [i~^q~ n+ ^ - ^q n+ ^, 
5 n = (n\8\n) = {w/z) fi^q n+ ^ - fi~^q~ n ~^, 

ii I ii — 1 n^i2n+2 ,, — 1„— 2n— 2 

In — \ n 7 \ n + 1) = U 2 q 2 : , n < N — 1, 

q-q- 1 

(3 n = (n\(3\n-l) = (w/z)(q-q- 1 )ti-*q- n+ v, n>0. (19) 
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Note that this representation can be extended to generic q but is then infinite-dimensional, 
i.e. n e Z> . The reason for this particular choice of the representation is explained in 
[HE]. Here we simply recall that the following exact sequence holds [l], 

- ^ K ® 7r?> ^ Tr^r 1 - 0, w = w'g" 1 = w"q = z/fi, (20) 

with the inclusion i and the projection r detailed in pQ (cf Section 4.1 and 4.2). From 
this decomposition of the tensor product one now derives a functional equation of the 
following type 

T{z)Q,(w)=<P l (z)Q m (w') + ct )2 (z)Q w -,(w'') (21) 

where 

N' 

T={Tr®l H )T = A + D and Q M = (Tr <g>l«)Q = V Q nn . (22) 

^ n=0 

Here <pi,4> 2 are some coefficient functions whose precise form depends on the choice 
of the quantum space 7i which is as yet unspecified. For the XXZ spin-chain we will 
present them below. 

2.2 Infinite-dimensional auxiliary space for q N ^ 1 

Rossi and Weston introduced in ^§] the following infinite-dimensional four-parameter 
representation n + = n + (w; s , si, s 2 ) of the upper Borel subalgebra U q (b + ), 

7r+(ei)|n) = |n-l), 7r + {q hl ) \n) = s g" 2n \n) , ir + (q ho ) = ir + {q~ hl ), 

vr + (e ) \n) = fsis 2 (g ~ g ^ + . W — + Sl q 2n + s 2 q- 2n ) \n + 1) . (23) 

Here nfZ and the representation space is thus infinite-dimensional. For convenience 
we rescale the parameters according to 

si,2 -> t-1,2 = -(<? - g _1 ) 2 w _1 s lj2 , s ^r = s . (24) 
The matrix elements for the intertwiner L are now calculated to be [TS] 

a n = (n|a|n) = (w/z) r 2 r 2 g~ n+2 — r 2 g™, 

i i 

5 n = (n\6\n) = (w/z)rirQq n — rQq~ n , 

_ i 

7 n = (n|7|n + l) = (? - ) r o 2 g" +1 , 

(3 n = (n\/3\n - 1) = (iw/z) r 2 g n+1 _ . (25) 

Similar to the root of unity case one has a decomposition of the tensor product 7r + ® tt^ 
according to the exact sequence [Hi] 

-> vr + (wg 2 ;r + ) ^ vr + (w;r) <g> vr^ A 7r + (wg~ 2 ;0 -> 0, iu = 2, (26) 



6 



with 

r = (r ,r 1 ,r 2 ) and = (r g ±1 , r x q T2 , r 2 ) . (27) 

The inclusion and projection map now are [Til] 

i\n) + =r - p— In) ® |0) + In - 1) <g> |1) and r In + 1) ® |0) = |n)~ . 

q-q 

Again this representation theoretic fact is the platform for deriving a functional equation 
of the type 

T{z)Q{z; r) = r + ) + V 2 (*)Q(*; r"), (28) 

where one now has 

oo 

Q(w; r) = ( Tr ®1 H )Q = V Q nn . (29) 

7r+(w;r) ' * 

n=— oo 

The reader will have noticed that the sum in the definition of the auxiliary matrix is 
now infinite. In fact, this can cause convergence problems. From a mathematical point 
of view one should therefore treat q as an abstract indeterminate (rather than a complex 
number) and view the matrix elements of ()29|) as formal power series in q. We will return 
to this point later when we calculate the spectrum of (J2HJ). See also the discussion in 

EH]- 

2.3 Commutation relations with the Yang-Baxter algebra 

Having specified explicitly the matrix elements of the intertwiner (fTTj) at and away 
from a root of unity we can exploit the Yang-Baxter equation (JT5j) in order to derive 
the commutation relations of the matrix elements (|14|) with the generators © of the 
Yang-Baxter algebra. One finds 

akQkiA = aiAQ k i + j l _ 1 BQ k i^ 1 -P k Q k _ u C, (30) 

a k Q k iB = 5iBQ H + (3 l+1 AQ k i +1 - (3 k Q k -uD, (31) 

SkQkiC = aiCQ k i + 'Y l _ 1 DQ k i- 1 -j k Q k+ iiA, (32) 

hQuD = 5 t DQ kl + (3 l+1 CQ ki+ i - <y k Q k+ll B . (33) 

Here we have suppressed the dependence on the spectral variables in the notation which 
for say the second identity is explicitly given by 

Q kl {w)B{z) = 

$i{w/z) (3 l+1 (w/z) /3 k {w/z) f \r>( \ 

— — - - B(z)Q k i(w) + A{z)Q kl+l (w) -— - Qk-xiiw) D(z) . 

a k {w/z) a k {w/z) a k {w/z) 

Depending on the type of representation determining the auxiliary space of Q there are 
different boundary conditions. For example, if q N = ±1 then 

P Q = a QoiB = 8i BQ i + p l+1 AQ m+l , 
P N i = =>- a k Q k N'-iB = 5jv'-i BQ kN i_i — fi k Qk-iN'-i D . 
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Away from a root of unity the auxiliary space is infinite-dimensional and there are no 
boundary conditions. 

Notice at this point the difference in the alternative construction procedures for 
auxiliary matrices. In order to obtain an auxiliary matrix which commutes with the 
transfer matrix we used the concept of intertwiners leading to p5j) . The latter implies the 
commutator [T(z),Q(w)] = which corresponds to the following non-trivial equation 
in terms of the Yang-Baxter algebra, 

f lk 1 £>Q kk ^ l _ 1±_ Q k+lk B] = ( — Qk-ikC — ^ k+l CQkk+i ) • (34) 
V V ak 6k ' k ^ ak 6k ' 

In contrast Baxter's method described e.g. in ^H] relies on the "pair propagation through 
a vertex" property to construct an auxiliary matrix which commutes with T . The 
relation ()15j) . which is the key to relating the Yang-Baxter algebra (JUJ) with the matrix 
elements ()14|) . is missing in Baxter's approach. 

3 Action on Bethe states 

Since the transfer matrix and the Q-operators commute, they must allow for a common 
set of eigenstates. This motivates us to compute the action of the Q-operator on the 
Bethe eigenstates of the transfer matrix, i.e. we now make the crucial assumption that 
the quantum space TL contains a vector |0) w , called the pseudo- vacuum, which satisfies 

A |0> w = |0> w (0\A\0) H , C |0> w = 0, D |0> w = |0) w (0|D|0> W , (35) 

and 

Qkk\0) n = \0) n (0\Q kk \0) H , Q jk \0) n = 0, 3 >k. (36) 

In fact, when the quantum space Ti carries a representation 7r w of the quantum group 
Ug(sl 2 ) or U q {bJ) this pseudo- vacuum is identified with a (unique) highest weight vector. 
The above properties then follow from the intertwining property © of the monodromy 
matrices. For instance, one has for Q that 

[Q, <{q h ) ® 7Tn(q h )} = =► ^ H {q h )Q 3 ^ H {q~ h ) = q 2 ^Q ]k , 

which then implies (}3l)j) by exploiting the fact that \0) n is highest weight. 
We now define a "proper" Bethe eigenstate as a vector of the form 

n B 

\z 1 ,...,z nB ) n = Y[B(z j )\0) H , (37) 
i=i 

where the parameters Zj = Zj(q, A) are finite solutions to the "generalized" Bethe ansatz 
equations 

(0\A(zi)\0) H q nB P B {z t q- 2 ) + (0|D(^)|0) W q- nB P B { Zl q 2 ) = . (38) 
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Here we have introduced for convenience the "Bethe polynomial" 



'ID 



Pb{z) = Y[{1 - z/z-) with Zj = Zj(q, A) . (39) 
i=i 

This set of coupled non-linear equations is sufficient to guarantee that the Bethe states 
()37p are indeed eigenstates of the transfer matrix T = A + D. The proof follows the 
same lines as the well-known computation for the XXZ spin-chain [20] . 



T(z) = (0\A(z)\0) n g Wfl P p ( 7"x 2) + <°P(*)|0>* V~ nB ^TT- ■ ( 4 °) 



We restrict ourselves to finite solutions of the Bethe ansatz equations (|3s|) in order 
to take into account certain peculiarities which can occur at roots of unity. Some of the 



finite Bethe roots Zj at q 7^ ±1 can tend to zero or infinity when the root of unity 
limit is taken. In this limit it might also occur that a subset of Bethe roots {<^}i e z , 
forms a complete string, 

lim J] (z - z lt {q')) =]!(*- = zN ' ~ Z M N '> I*' = ±i> ( 41 ) 



and so drops out of the equation (|3q) [22] . In terms of the Yang-Baxter algebra the 
occurence of a complete string corresponds to the vanishing of the following operator 
product, 

lim J] B(z u tf)) = 0, q N ' = ±l. (42) 

q'— *q 

Hence, the Bethe states (j3Tj) evaluated at a root of unity might not yield a complete set 
of eigenstates. Moreover, the fact that the transfer matrix and the Q-operator possess 
a common set of eigenvectors does not necessarily imply that the Bethe states are also 
eigenstates of the Q-operator when degeneracies are present. This is precisely the case 
when q N = ±1. For instance, in [T|, auxiliary matrices have been constructed for the six- 
vertex model at roots of unity which do not preserve the total spin and whose eigenstates 
therefore are different from the Bethe states. We will return to this discussion when 
we specialize our general setup to the XXZ spin-chain. For the moment we keep the 
calculation as general as possible. 

Let us start with the simplest case, only one "Bethe root" is present. That is, we 
evaluate the B operator at solutions zo to the equation 

(0|i4(2b)|0) w =<0|£>(2b)|0) w • (43) 

In order to compute the action of the Q-operator on the corresponding eigenstate of the 
transfer matrix we use the relation 



9 



which follows from the identities (|3U|h (|3ip. Note that if = / all of the rightmost 
operators on the right-hand-side of the above equation possess the pseudo-vacuum as 
eigenvector. Taking the trace in fi4*|) we obtain (the boundary terms work out correctly 
ifg*' = ±l) 

QkkB = 

k 

bJ2(-- Qkk+y^— Q k -i k (A -d)+j2 q*-i k+ i c . 

^ \a k a k a k+1 J k ak k ak ° k + 1 

When acting on the pseudovacuum the second and third term on the RHS of the equation 
vanish, which leaves us with the eigenvalue 

Q(w)B(z ) \0} H = \j2 (— ~ (0\Qkk\0} n \ B(z ) \0) H . (45) 

Here the sum runs over Z^/ when q N = ±1 and over Z when q is not a root of unity. The 
calculation for Bethe states with multiple Bethe roots follows the same logic, though it 
is now much more involved to show that all of the "unwanted" terms vanish due to the 
Bethe ansatz equations; see the appendix. 



3.1 Conjecture 

For a Bethe state with % Bethe roots we conjecture the following formula for the 
eigenvalues of the respective auxiliary matrices: 

Q(w)Y[B( Zj )\0) n = (46) 

(e <oiq»Min> w n ( ~4~t4 - Vr/f )7 ffl ) I ii B fe)io>« 

(V fJiKMw/Zj) OL k (w/zj)a k+1 {w/zj)J J A± 

Note that the particular choice of the quantum space only enters via the pseudo- vacuum 
expectation value (0|Qfcfc|0) w and the Bethe roots {ZjYj=\- The combination of matrix 
elements appearing in the product of the eigenvalue expression only depends on the 
auxiliary space. The above conjecture can be verified for states with % = 2,3 by 
employing the intermediate steps detailed in the appendix. For n B > 3 we will provide 
further support by showing that the eigenvalues satisfy the functional equations (}2"Tj) 
and (|28|) respectively. For this purpose we need to determine the coefficient functions in 
f)2ip and (|28|) first. The latter depend on the choice of the quantum space and we now 
specialize to the XXZ spin-chain. 
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4 The XXZ spin-chain 



As remarked upon earlier, one chooses for the familiar case of the inhomogeneous six- 
vertex model or XXZ spin-chain with quasi-periodic boundary conditions, 

M 

^ = (gK™ and n = (c 2 f M - (47) 

m=l 

Here {C m } are some complex inhomogeneity parameters. The transfer and auxiliary 
matrix can then be written as the trace of an operator product 

T(z) = Tr X^Rvm^/Cm) ■ ■■Roiiz/Q, (48) 

and 

Q(w) = Ti\<^ l L 0M {w/C M ) ■ --Loiiw/Q ■ (49) 

Here R and L are the i?-matrix and the L-operator specified earlier. Notice again the 
difference with Baxter's construction procedure [T^]. His final Q, which commutes with 
the transfer matrix, is in general not of this simple form. For the pseudo-vacuum and 
the associated "expectation values" one finds 

M . 

|0>„ = |0) c2 ® • • • ® |0) c2 , (0|A|0> W = A, (0|£>|0> w = \-\ M H ^TZ7- . ( 50 ) 

m=l Zq ^ m 

and 

M M 

(0|QH|0> w = \<W 91 n <*(w/U =► Q nn \0) n = \- 2n J] a n (w/U \0) H . (51) 

m=l m=l 

Here we have slightly modified our earlier conventions. Instead of taking the Cartan 
generator h in the exponent of the twist parameter, we introduced for convenience hi 
which in the two representations (J 18)) and (|23|) is given by 

7r£(/i') |n) = -2n \n) and n + (ti) \n) = -In \n) . (52) 

This simply amounts to a renormalization of the respective auxiliary matrices by an 
overall factor. We continue to treat the two cases of q being a root of unity and q not 
being a root of unity separately. 

4.1 Roots of unity 

When q N = ±1, the auxiliary space is finite and upon inserting the expressions (|T§j) 
into the conjectures formula (jH|) one obtains 

M 

(0\Qnn(w)\0) n = \- 2n »f q Mn+ f n(^~V 2 7Cn " 1), (53) 



m=l 
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and 

n B 

Q,(w)Y[B( Zj )\0) H = (54) 

2=1 

Here we have used the relation between total spin and the number of Bethe roots, 
2S Z = M — 2hb- We emphasize that the Bethe roots Zj are assumed to be finite 
solutions to the Bethe ansatz equations at a root of unity. This is different from the 
polynomial which is obtained by solving the Bethe ansatz equations away from a root 
of unity and then taking the root of unity limit. As pointed out before, in this limit one 
might encounter vanishing or infinite Bethe roots as well as complete strings. Before 
we investigate these issues we first check the eigenvalues of the auxiliary matrix for 
consistency by making contact with the functional equation (|21|) . 



4.1.1 The TQ-equation 

Using the explicit form of the inclusion and projection map one calculates for the twisted 
XXZ spin-chain the following coefficient functions in (|21|). 

T(z)Q lx (z/fi) = A" 1 ?" ( f[ V_ C ; ) QMn) + M f Q M -i (zq-'/fi) . (55) 

\rn=l Zq L ' m J 

Employing the known expression (|40|) for the eigenvalues of the transfer matrix within 
the framework of the algebraic Bethe ansatz [20J, 

one verifies that the conjectured eigenvalues of the auxiliary matrix satisfy (|55|). Notice 
that we have implictly made the assumption that all Q-matrices in (|55|) commute with 
each other. For A = 1 this has been proven in PQ. Using the explicit form of the inter- 
twiners employed in this proof one verifies that the same holds true for quasi-periodic 
boundary conditions as long as A w = 1. This is sufficient to cover the symmetries 
investigated in 0. 

4.1.2 The degeneracies at roots of unity q N> = ±1 

Recall from |S] that due to the loop symmetry at a root of unity q N = ±1 the 
eigenspaces of the transfer matrix are organized into multiplets containing states whose 
total spin S z varies by multiples of N'. The Bethe states (J37|) with finite Bethe roots are 
assumed to correspond to the highest weight states in such multiplets, i.e. acting with 
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the generators of the respective symmetry algebra on this state one successively obtains 
the whole degenerate eigenspace of the transfer matrix. (See |23*1 |2U 12] for examples of 
the highest weight property) The states within a degenerate multiplet are not of the 
form (jHZj), whence our result (JH4"J) only applies to the highest weight vectors containing 
finite Bethe roots. Therefore, it does not yield the complete spectrum of the auxiliary 
matrices Q^. We now compare for these highest weight states the identity (jSljl with 
the formulae (11) and (19) in [2] which have been proven for N = 3 using functional 
relations and conjectured to hold true for N > 3 based on numerical data. 

Recall from 2J that the eigenvalues of the auxiliary matrices for the homogeneous 
chain, {( m = l}£f =1 , and periodic boundary conditions, A = 1, were shown to be of the 
following general form (see formulae (11) and (19) in j2j), 

Q,(w = z/fi) = N ll z^P B {z)P B {z^ 2 )P s {z N ' ,fi 2N ') . (57) 

In order to match our result ()54|) from the algebraic Bethe ansatz computation with the 
formulae (11) and (19) in [2] we have to identify 

ns 2kS z I —2k i \M 

u^p s (z"',^ i) = nd - = e £ ZjLzX-i, ■ <«) 

3=1 keZ N / 

Note that in our derivation of (|54jl we have assumed infinite and vanishing Bethe roots 
to be absent. Thus, we have to set = in formula (19) of j2J. Numerical evidence 
suggests that this is only true in the commensurate sectors 2S Z = OmodiV. One then 
finds for N = 3 agreement between the results in [2] and our present calculation up 
to a trivial redefinition of the normalization constant M^. Recall from the discussion 
in [2] that the rational function ([58)1 is in fact a polynomial* in the variable z N whose 
roots a,j contain all the essential information on the irreducible representation of the 
loop algebra spanning the degenerate eigenspaces of the transfer matrix. The reader is 
referred to [2] for details. 

In contrast to the transfer matrix the auxiliary matrices stay non-degenerate and 
in addition to the eigenvalues corresponding to the highest weight state one also needs 
to compute the remaining eigenvalues of the auxiliary matrix within the multiplet. The 
examples explicitly worked out for N = 3 and M = 5, 6, 8 in |2j suggest the following 
picture: 

Set A = 1, {( m = l}m=i an d assume we have a multiplet in a commensurate sector, 
i.e. the highest weight state has spin 2S Z = OmodiV. Then the form of the eigenvalue 
(p>4"|) remains the same except for a change in the //-dependence of the polynomial ([58)1 . 
The examples in [2] showed that this polynomial changes within the multiplet according 

*We briefly recall the argument. First, notice that the Bethe ansatz equations 1)38(1 ensure that the 
rational function l(58|l is pole free in z. Hence, it must be a polynomial in z. Secondly, the function 
B58J1 is obviously invariant under the replacement z — + zq 2 . Thus, any zero z' ^ occurs inside a string 
{z'q 2e }g£z N , and gives rise to a factor (z N — z' N ) in the root decomposition of the polynomial. 
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to the following rule, 

S z : P s {z N ',fi 2N ') = Y[{l-z N 'fi- 2N '/a J ) (highest weight state) 

3=1 

n s 

S z -N> : P s (z N \ = (1 - z"M f[(l - z N '^ N '/a 3 ) 



3=2 

S* - 2N> : P s (z N '^ 2N ') = (1 - z N '/ ai )(l - z N '/a 2 ) - z N ' ^ 2N ' / a,] 

J'=3 



-S 2 : P S (^') = IJ^ - zN '/ a j) ( lowest wei S ht state ) 

3=1 

That is, as one steps through the multiplet the polynomial (J58j) successively looses 
factors of fi~ 2N until one reaches the lowest weight state, where it does not contain any 
/i-factor. We can confirm this picture by computing the eigenvalue corresponding to 
the lowest weight state. Invoking the transformation of the auxiliary matrix under spin 
reversal 9^ = [] m=1 a x m one finds (cf formula (39) in 0), 



n B n B 

n(v.\ ln\ 

H 



Q^w = z/n)Y[C(z J )m\0} H = const, x *RQ M -i(w = z/fi) JJ B( Zj ) |0) 

3=1 3=1 



In accordance with the picture outlined before, the sum yielding the polynomial Ps is 
now independent of the parameter /i. 

4.1.3 Infinite and vanishing Bethe roots 

Besides the missing states within a degenerate multiplet of the transfer matrix, we are 
also missing those states which contain infinite or vanishing Bethe roots in the root of 
unity limit. A concrete example for the homogeneous chain, i.e. {( m = l}^f =1 , is given 
by setting M = 5, S z = 1/2 and considering the subspace of momentum P = 0. There 
one finds for 2A — q' + q'~ x with q' not a root of unity the following Bethe roots 



h = b? 1 = ^(l + A - a/5 + A(A-2) + \j (l + A - a/5 + A(A - 2] 



2 - 16 



Here we have temporarily introduced the parameterization bi = q'(l — Zi)/(1 — Ziq' 2 ) in 
order to accommodate the occurrence of infinite roots. Taking the root of unity limit 
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q' — > q = e 27rJ / 3 one now computes 



lim T(z) = 1 + ^] , lim b x = lim b 2 1 = e 2 ™ /3 = q, Zl = 0, 2 2 = oo . 

Therefore the /mite solutions to the Bethe ansatz equations do not give a complete 
set of eigenstates at roots of unity. Notice that it might even happen that the highest 
weight state in a multiplet contains infinite or vanishing Bethe roots. 

4.2 The case when q is not a root of unity 

A priori one might expect that things are now simpler than the root of unity case as the 
only degeneracies which occur in the spectrum of the transfer matrix are due to spin- 
reversal symmetry. However, the auxiliary space is now infinite and one has to deal with 
potential convergence problems when taking the trace in (J29|) . This becomes evident if 
we insert the expressions (12 5|) into the conjectured formula for the eigenvalue expression 
of the auxiliary matrix. One obtains for the pseudo-vacuum expectation value 

M M 

(0\Qnn(w)\0) H = ^ II «>/Cm) = A^r' V f U (wr 2 q- 2n+2 / ( m - 1) (59) 

m=l m=l 

and the action of the auxiliary matrix on a Bethe state produces 

n B 

Q(z; r , n, r 2 ) JJ B(zj) \Q) n = 
i=i 

{\-2£ 21S Z y\ m ( ?r n - 2i + 2 1 Y ^^^ UB 
r^P B ( Z )P B ( Z r ir2 ) X l ( Um ±l^ ( ,C Z: II B &) \ Q )n ■ 
PB[zr 2 q M+2 )PB{zr 2 q M ) J jj^ 

The above expression for the eigenvalue might not be convergent. In fact, it was dis- 
cussed in [Hi] for 7*1 = r 2 = and A = 1 that the matrix elements of the Q-operator 
contain the formal power series 

S( q 2 n = J2q 2lS \ S* = f-n B (60) 

in the deformation parameter q. If one removes this factor in the sector S z = by an 
ad-hoc renormalisation the remaining auxiliary matrix can be identified with Baxter's 
expression (101) in jTHj. See the discussion in Section 5 of [19J for details. Obviously, 
our algebraic Bethe ansatz analysis of the auxiliary matrices reproduces these findings. 
In the limit r\,r 2 — > we obtain apart from the factor (|HU|) and the constant r^ S2 , the 
Bethe polynomial ()39|) as eigenvalue as we should. However, the required renormalisation 
is certainly an unsatisfactory answer to the convergence problem, in particular as the 
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factor (jHOJ) is needed outside the spin-sector S z = in order to satisfy the functional 
equation (|28|) which for the XXZ spin-chain is calculated to 



Q(z;r ,r 1 ,r 2 )T(z) = (61) 



M 



A 1 q^Q{zq 2 ;r q \ r x q 2 , r 2 ) + \q * ( TT * -f ) Q{zq £ i ;r q,r 1 q 2 ,r 2 ) . 

In order to obtain a convergent expression let us modify the auxiliary space by restricting 
it to an invariant subspace at special values of the parameters ri )2 . 

4.2.1 Truncation of the auxiliary space 

In order to establish convergence it would be helpful to restrict the sum in the eigenvalue 
expression of the auxiliary matrix either to the positive or negative integers and then 
choose the twist parameter appropriately. In order to achieve this we have to find values 
for the free parameters (r , r%, r 2 ) such that the representation (|23|) contains an invariant 
subspace. To this end we need for some m e Z the truncation condition 

rir 2 + 1 — nq 2m ° — r 2 g" 2m ° 
e Imo}^ = w _ 2 \m )^ + = . (62) 

This is easily achieved by setting 

either n = q~ 2m ° or r 2 = q 2m ° . (63) 

For simplicity let us choose m = and r 2 = 1. Denote by 7r< = n<(w ; r , ri) the 
irreducible subrepresentation of (|23|) which is obtained by restriction to the negative 
integers (including zero) and set 

o 

Q<(w; r , n) = (Tr ®1 W )Q = ^ Q n „ . (64) 



n=— oo 



Since the parameter r 2 does not shift in the decomposition ()26|) this truncation is con- 
sistent with the functional equation (J2*8j) . albeit we have to supplement the projection 
map r by defining 

T \ ),i( W ;r , ri ) ® = g2 _ x l°>^(™rVo«rW) r x ^ 1, w = 2 . (65) 

For the semi-infinite representation space we now obtain the following eigenvalue asso- 
ciated with a Bethe state (JHZJ), 

ostein) = x-wp^p^z) p a v- g^ff^V^ . fee) 

Here we have denoted the eigenvalue by the same symbol as the operator. Unlike in the 
previous expression the sum now only ranges over the positive integers. This puts us 
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in the position to prove absolute convergence of the series in (J6bJ). We treat the cases 
|g| = 1 and |g| ±:L > 1 separately. 

Suppose q is of modulus one, then we deduce from 



1= E 



\2l-MS* Ylm=l( Z( l 2e /Cm 1, 

A g 



P B (zq*)P B (zqV-*) 



- ipwy 11 <67) 



that absolute convergence is guaranteed as long as |A| < 1. 

Now let |g| ±:L > 1 and set q = expt>, z = exp2u, ( m = exp2£ m . Rewriting 



oo 

E 

1=1 



\2t-2lS z Ylm=l( Z( l 2e / Cm 1) 

A q 



P B (zq 2i )P B (zq 2i - 2 ) 

*\"\ m' lB \ Z \~~ nB ]lZ=l l Slllil l»-<;mT^IIWI ; < congt I x |2<| |±2<|S* 

n ni3 I ainU« - 11, -U A,"! sinlif7/ - 11, -U f/> — _ ' 1 1 |y| 



2|A 


2f 


Is 


n B 


Z 


f-n B U M 

L Lm=l 


sinh(u-£ m + fo)||C m 


i 

2 




1 


sinh(w 


— Uj + tv) sinh(-u — Uj + (£ — l)v)\/ 


\ Z j\ 



we find absolute convergence for |A| < \q\ ThI ^ 2 . Thus, in summary we are left with the 
condition 

|A| < \q\™'\ \q\ ±l > 1 (68) 

which now includes the case when q is of modulus one. 

Having assured convergence we can now verify that the eigenvalues (J66|) satisfy the 
functional equation (J61|) . Doing so we again implicitly made the assumption that the Q- 
operators in the functional equation commute with each other. Since the auxiliary space 
is still semi-infinite, the proof of commutation via the construction of the corresponding 
intertwiners (similar to the root of unity case) is less feasible. Instead we are going to 
exploit the completeness of the Bethe ansatz when q is not a root of unity. 

Provided one accepts that the Bethe states (|37|) and their counterparts under spin- 
reversal 

n B 

^\z 1 ,...,z nB ) n ^Y[C(z j )^\0) n , «H|0) W = |l) ca ® |l) ca , (69) 

3=1 

form a complete set of eigenstates which span the whole quantum space TC, one deduces 
that the three-parameter family {Q<(z; vq, n)} of auxiliary matrices can be simultane- 
ously diagonalized as the Bethe roots only depend on q and A. Hence, we must have 
[<2<(2;ro,ri),<2<(w;r(,,ri)] = 0. 

Notice that the commutation of the Q-operators for arbitrary spectral parameters 
together with the explicit form of the matrix elements (J2*5j) implies that the eigenvalues 
(JBTIj) are polynomials of degree < M in the spectral variable z. Unlike the root-of-unity 
case, however, the sum by itself is not a polynomial. One now has to include a factor 
P B in front of the sum in order to cancel the simple poles of the denominator at z — Zj. 
Employing the Bethe ansatz equations (JH8|) one then deduces that for any closed contour 
Ct around Zjq~~ 2e with integer £ > one has 

<l rl-y P (-A \'^ n -2eS z Ilm=l( Z 1 2e /Cm ~ ]) _ n 
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The vanishing of the contour integral signals that the above rational function is indeed 
a polynomial in z. This is consistent with (|66jh where the parameter t\ in the argument 
of the first factor is arbitrary. 

For completeness we briefly discuss the transformation under spin-reversal in order 
to obtain the eigenvalues of the auxiliary matrix "beyond the equator", i.e. for the 
states (|69|). Acting with the spin-reversal operator on the auxiliary matrix we obtain 

XQ<(z; r„, n)m = Tr \*±W )9l o* M L 0M (z/( M )* x M ■ ■ ■ cx^oi^/CiK . 

which amounts in terms of the intertwiner L to the replacement 

{a, /3, 7, 7, A a}. (70) 

From (|25|) one derives the following identities for the matrix elements of the intertwiner 
L, 

a n (wq~ 2 ,q;r ,r 1 ,r 2 ) = 5 n (w, q' 1 ; 1 , r 2 , n) 
5 n (wq~ 2 ,q;r Q ,r 1 ,r 2 ) = a n (w,q' 1 ;rQ 1 ,r 2 ,r 1 ) 

Pn1n-l\(wq- 2 ,q;ro,ri,r2) ~ Pnln-l I (w ,q- 1 ;r^ 1 ,r 2 ,r{) V"U 

Notice that only the above combinations of matrix elements contribute to the trace of 
the monodromy matrix ([13)1 . i.e. the auxiliary matrix (|49p. Hence, we have the identity 

^Q<(z<r 2 , g ;r ,r 1; r 2 = 1)01 = Q<(z, q' 1 ; r^\r 2 = l,n) . (72) 

Thus, spin-reversal leads to the consideration of the truncated auxiliary matrix with n = 
1 and r 2 arbitrary. We now discuss this case in the context of g-oscillator representations. 
Before we have the following 

Remark. In order to ensure convergence for |g| = 1 we have excluded the quasi- 
periodic boundary conditions with |A| > 1. However, in order to recover for example 
periodic boundary conditions, A = 1, one may proceed as follows. Employing (|61j) 
one can express the transfer matrix eigenvalues in terms of the auxiliary matrix at 
|A| < 1. While at the moment we do not have an analytic argument due to the implicit 
dependence of the Bethe roots on A, numerical calculations suggest that one can then 
safely take the limit A — > 1 in order to recover the eigenvalues of the six-vertex transfer 
matrix with periodic boundary conditions. 

4.2.2 g-oscillator representations 

We now simplify the representation (|23J) of JH] further in order to make contact with 
the results in jT7] for the Coulomb gas formalism of conformal field theory. We explicitly 
show that one can derive the analogue of the functional equations therein also for the 
XXZ spin-chain. 

The g-oscillator algebra is defined in terms of the generators 

qe + e_ — q~ 1 e_e + = (q — q^ 1 )^ 1 and [h,e±]=±2e±. (73) 
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We define two representations q ± by identifying eo — > e±, e\ — > e T , hi = —h ^ ±h 
and setting 

£ + = 7r + (w; r = 1, n = 1, r 2 = 0) and g_ = tt + (w; r = 1, r x = 0, r 2 = 1) . (74) 

Again we notice that the representation space in truncates and it will be understood 
that the representation spaces associated with g± are the modules obtained by acting 
freely with e± on the highest (lowest) weight vector |0). Let Q ± = (Tr e ± ®1)Q then we 
find for the eigenvalues associated with a Bethe vector (JHZJ, 

<t{z) = {-l) M P B {z) £ \ 2i q~ 2es * = , (75) 



and 



Q (z) - P fl (z) ^ A g p B( ^ +2)Pfl( ^ ) • (76) 



£=0 



For absolute convergence we have assumed that ()68|) holds as befrore. Notice that in 
order to satisfy the functional equation (}6*T]) with the transfer matrix we have to keep in 
mind that the operators Q still depend implicitly on the parameters (ro,ri,r 2 ) which 
shift in (J26j) . This is particularly important in the case of Q + = Q<(r = l,rj = 
l,r 2 = 0) where in the functional equation the shift T\ — > riq ±2 occurs and, hence, the 
truncation condition (Jo^j) for the auxiliary space changes. 

This highlights the importance of introducing the free parameters (tq, r\, r 2 ) in (|23|). 
Even if one would directly start the discussion with the simpler looking representations 
(jUj) the decomposition of the tensor product (J2EJ, which underlies the TQ-equation 
(|28p. leads to the consideration of the more general representation ()23|). (See also the 
comments in the introduction of [lj for the root of unity case.) We should therefore 
understand the expressions (|75j) and (fTUj) as a decomposition of the general solution 
(|66p. Namely, we can write (JBTIj) as a product of ([75]) and ()76|) (compare with formula 
(4.10) in [T7j), 

Q<(s;r ,ri) = (-l) M r - 5Z (l - X 2 q- 2SZ )Q + {zr l )Q~{z) . (77) 

One now easily derives the analogue of the functional relations reported in |T7j for the 
Coulomb gas formalism. For instance, we find for the eigenvalues the following formula 
corresponding to the "quantum Wronskian condition" (cf equation (4.3) in [T7|). 



en :<r)Q-(--) aV 2 - s '-Q + (--)Q-(^ 2 ) = U ™: l( \ 2 ^{ Cm) ■ <7> 

1 — A g - ^ 



This is a special case of a more general relation (see equation (189)) below) involving the 
fusion matrices of the six- vertex model, which we discuss next. 
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5 The six-vertex fusion hierarchy 

In order to provide further support for our algebraic Bethe ansatz computation of the 
eigenvalues of the auxiliary matrices, we now make contact with what in the literature 
is known as the fusion hierarchy. We follow a similar line of argument to that in for 
the eight-vertex model. First, we will briefly review the representation theoretic aspects 
in the construction of the fusion hierarchy which will be the key to connecting fusion 
and auxiliary matrices at roots of unity. We specialize at once to the XXZ spin-chain 

flUD- 

We start by introducing the monodromy matrices associated with the fusion matrices. 
Denote by ^-f 1 : U q (sl2) — > EndC n+1 the spin n/2 evaluation representation of the 
quantum loop algebra, i.e. 

7rW(ei)|m> = [n - m + l] q \m - 1) , 7r< B >(/ ) = z' 1 ^\e 1 ), 

^\fi)\m) = [m + l]Jm+l>, 7rf(e ) = «r< B) (/i), 

ttJV 1 )!™) = 1 n ~ 2m H, 4 n V )=^ n) (<T /ll ) ) (79) 
with m = 0, 1, n. Let 

X (n+i)( zg -n-i) = ®ir H )\ hm ReEndC n+1 ®H, (80) 

with 

Z> +1 ) ( w ) = U$ ® vrfii) R G End C n+1 ® C 2 . (81) 

Note that we have labelled the fusion matrices by the dimension of their auxiliary space 
instead of the spin. The matrices (JTj), (jSUj) and (}8~Tj) again satisfy the Yang-Baxter 
relation. The matrix elements of (j8Tj) w.r.t. the representation 7r[.=i are explicitly given 
by 

(0|L (n+1) (uO|0) = P + vr (rt) (^) -p_vr (n) (g-t), 

<0|L(" +1 )H|1> = p + (g-g-V n) (?V B) (/i), 

(l\L^ +1 \w)\0) = p^(q-q- 1 )7r^(e 1 )^ n \q-^), 

(l\L^ n+1 \w)\l) = p + vrW(g-t) _p_ 7r W(gf) (82) 

with the coefficients p ± satisfying the constraint p + / p_ = wq. From the following non- 
split exact sequence describing the decomposition of evaluation representations of the 
quantum loop algebra f/ g (s/ 2 ) |2*T] . 

-> 7r£ n r 1} A ttW ® ttW A tt^ -> 0, w = ^V 1 = <A = z<f +1 (83) 

one derives a functional relation, known as the fusion hierarchy (see e.g. [2ZII2E1EII and 
also 30J for an alternative form of the fusion equation). 

M M 

T^\z)T^\zq- 2 ) = T^ +1 \zq- 2 ) J] (zq 2 /C m - l) +T^ 1 \zq 2 ) J] (z/( m - 1) . (84) 

m=l m=l 
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Here 

T^ n+l \zq- n ~ 1 ) = ( Tr ®l)T (n+1) (^ n - 1 ) G End ft . (85) 
The coefficients in (I82j) have been chosen as 

P+ = wq, p_ = 1, (86) 

and we have identified 

M M 

TW(zq- 2 ) = q-*T(z) ]J (zq 2 /( m - 1), T®(z) = ]J (zq 2 /( m - 1) . (87) 

m=l m=l 



The coefficient functions in f!84|) have been calculated using the explicit form of the 
inclusion and projection map in 



7r (n x) 3 \m) ^ |m)' = [n - m] |m) <g> |1) - q n ~ m [m +1] |m + 1) ® |0) , 

and 

[n — m + 1J 

Inserting the explicit expression for the eigenvalues of the six-vertex transfer matrix 
from the algebraic Bethe ansatz (|4U|) we obtain the following formula for the eigenvalues 
of the fusion matrices associated with the Bethe states (|37jl . 

T<»> W = X—\^P B (z)P B (zi") ± AV^ g"^^,, 1 ) ■ (88) 

Here we have again denoted eigenvalues and operators by the same symbol. Pb denotes 
the Bethe polynomial defined in (|39J) whose zeroes are solutions to the Bethe ansatz 
equations (|38p. The proof follows via induction and is straightforward. Notice that the 
eigenvalues ()88|) are polynomials in z. The rational function in (|88p only has poles at 
2n which are cancelled by the factor PB{z)PB{zq 2n ) in front of the sum. The 
residue of the remaining poles is zero due to the Bethe ansatz equations. 

Away from a root of unity we can already connect the fusion hierarchy with the 
auxiliary matrices (f?3|) and ([76)1 . A straightforward calculation proves for the eigenvalues 
the relation 

C_i ) M \- 2 n 2SZ 

\- n q nS *Q + (zq 2n )Q~(z) - X n q- nSZ Q + (z)Q-(zq 2n ) = [ } Q T^\z) , (89) 

A f - \q~ b 

which corresponds to formula (4.1) in jTTj. The above functional equation not only 
holds for the eigenvalues but also for the operators provided the Bethe states constitute 
a complete basis of the quantum space. 
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5.1 Degeneracies at roots of unity 

The discussion of the fusion matrices has so far applied to the case of "generic" q (i.e. 
the deformation parameter being either a root of unity or not) and finite solutions to 
the Bethe ansatz equations entering the eigenvalues via the Bethe polynomial . We 
now specialize to the case q N = ±1 and start by showing that the fusion matrices 
exhibit the same infinite-dimensional non-abelian symmetries as the transfer matrix. 

Recall from [3J |S] that the symmetry generators for the XXZ spin-chain are given by 
(2 = 0,1) 

M 

Ef' ] = Bm^^AWCeD/W, 00) 

m=l 
M 

F t m = lim(g) 7 r«AW( / r)/[iV , ] 9 '!- (91) 

m=l 

Here A( M ) = (1 ® A)A( M ~V denotes the M-fold coproduct with A (2) = A andy is 
some number with q' N ' ^ ±1. The above form of the symmetry generators of [/(s^) is 
restricted to the commensurate spin-sectors 2S Z = mod N when A = 1 [3] but extends 
for the subalgebras U(b^) to all spin-sectors when A = q ±s * The extension of the 
loop symmetry from the fusion degree n = 2 to arbitrary n now simply follows by 
induction from the recursion relation ([84)1 . Thus, the eigenspaces of the fusion matrices 
also organize into multiplets containing states whose total spin differs by multiples of 
N'. 

Remark. In a similar manner one proves the loop symmetry of the higher-spin 
six-vertex models. The only major difference lies in the choice for the representation 7r^ 
defining the quantum space which changes to h-h = (5§™_i ■ An explicit calculation 

_____ ^-"^ III L ^ m 

shows that the analogue of the operators (|90j) is well-defined. Adopting the proof in 
jS] for the loop symmetry of the "fundamental" fusion matrix one shows again by 
induction that T^ n+1 ^ exhibits the same symmetries. This provides an alternative proof 
to the one given in |2*3] . 

5.2 Spectrum of the fusion matrices at q = 1 

In order to elucidate the spectrum of the fusion matrices at roots of unity let us discuss 
the limit q' — >• q from transcendental or irrational q' to q being a root of unity. In order 
to simplify the notation we will often denote this limit by the symbol q N — > 1 in the 
following. In order to perform the root of unity limit we need the explicit dependence 
of the Bethe roots Zi on the arbitrary deformation parameter q' . Unfortunately, this 
dependence is in general not known. However, for small chains up to the size M < 8 one 
can compute some Bethe roots analytically. According to these examples the following 
scenarios might occur in the root of unity limit: 

1. A Bethe root tends to zero, 

lim Zi = 0. (92) 
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We shall denote the number of such Bethe roots by no- 

2. A Bethe root tends to infinity, 

lim Zi = oo . (93) 

We shall denote the number of infinite roots by n^. 

3. There are N' Bethe roots which form a complete string 

!i m (^o) z hi • • • i z i N i_i) = i z ioi z h = z ioQ ■>•••■> Z i N i_ 1 = ■ z «o ( ? ) • (94) 
q" — *1 

Note that these complete strings obtained in the root of unity limit do not coincide 
with the zeroes of the classical Drinfeld polynomial (p)B*j) contained in the spectrum 
of the auxiliary matrices. 

Because of these three possibilities we deduce that taking the root of unity limit of 
the Bethe polynomial (|39|) is not the same as the Bethe polynomial containing the finite 
solutions of the Bethe ansatz equations at a root of unity, i.e. 

lim P B (z, q') ^ P B (z, q) with q N> = ±1 . 

In the presence of vanishing Bethe roots the root of unity limit of the Bethe polynomial 
might not even be well-defined. However, the eigenvalues of the fusion and transfer 
matrices only contain ratios of these polynomials. From the simple identities 

lim , ( ! ,~ Z ^ q2 l~ Z%) , = g"+ 2 , (95) 



Zi 



+o (zq 2i — Zi)(zq 



and 

lim (^ ~ z i)i z (f n ~ z i) = l (96) 
z z ^co (zq 2E — Zi)(zq 2e ~ 2 — Zi) 

we infer that the spectrum of the fusion matrices at a root unity changes to 

kmT (*)-A q P B (z)P B {zq )^A q • ( 97 ) 
Here we have set 

s = 2n + S z modN' . (98) 

Again we remind the reader that P B is now the "reduced" Bethe polynomial, i.e. it only 
contains the finite Bethe roots at q — 1 and no complete strings. From the identity 
(19 7j) it follows by a direct calculation that one has the following simplification of the 
fusion hierarchy in terms of the eigenvalues 

M 

lim T^ N ' +l \z) = (\- N 'q N ' s + X N 'q- N ' s ) T](zq 2 /( m - 1) + lim T^W) . (99) 
This formula has been reported for A = 1, ( m = 1 before in the literature |3T } I29 | IT^j. 
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5.3 Connection with auxiliary matrices at roots of 1 

We will now argue that the fusion matrix of degree N' at q N = ±1 can be identified 
with a special limit of the auxiliary matrix (|49jl. First, we observe that the corre- 
sponding auxiliary spaces have the same dimension. From the representation theory 
of quantum groups at roots of unity EB] one now deduces the following. Restrict 
the respective evaluation representations (fTSj) and (J79)l of the quantum loop algebra 
to the finite subalgebra U q {sl<2). If the values of the central elements in the respective 
^(•^2 ^representations determining the auxiliary spaces are equal then they must be 
isomorphic. The values of the central elements at a root of unity are 

T m : 7r^(ef')=7rW(/f) = 0, ^ N '\q h f' = ^ N '~ x) 

Q, ■■ ^(ef) = ^Uf ) = 0, ^(q h ) N ' = q- N 'li- N ' 

and for the Casimir c = q h+1 + q~ h ~ 1 + (q — q~ 1 ) 2 fe one finds 

TW : tt^(c) = q N ' + q~ N ', 

Qli ■ 7T M (C) = /i + Zi" 1 . 

Hence, in the limit /x — ► q N = ±1 both representation become isomorphic and one can 
therefore construct a gauge transformation in the auxiliary space rendering the fusion 
and auxiliary matrix equal (see [I] for details). Thus, we conclude 

lim T {N '\z) = lim Q„{z/fi) . (100) 



A" 



Note that in this limit the auxiliary matrix therefore also becomes degenerate. Never- 
theless, we can check the eigenvalues for consistency be comparing the expressions for 
the eigenvalues derived from the algebraic Bethe ansatz with the ones from the fusion 
hierarchy. One finds for the eigenvalues of the fusion matrix 

lim T^'\z) = \^'-^>P B ( Z )P B (zr') £ X2eq ^^f{tr 1] • 

^ t P B [zq M )P B (zq M 2 ) 

(101) 

In the commensurate sectors 2S Z = mod N where uq = = we can set s = in 
(jlOlj) . Then the spectrum of the auxiliary matrices derived for finite Bethe states at 
a root of unity coincides for fi = q N = ±1 with the one of the fusion matrix. Outside 
the commensurate sectors it has been argued (and proven for N = 3) in 2\ that the 
auxiliary matrices (|49|) have eigenvalues of the form (|57|) . When comparing with formula 
(11) in [2] one has to make the replacement — > in the notation. Here Ps denotes 
the previously discussed polynomial (}58|) whose zeroes are fixed by the representation 
theory of the loop algebra up to an factor /i 2JV , 

AT' 1 . if n x p / N' = N\ = (N'+l)s V- <£^ fL=l ( WCm ~ 1) (m) 



24 



The power of the monomial must be an integer by construction of the auxiliary matrix 
and the normalization chosen in ([19)1 . Interestingly, numerical examples show that 
vanishing and infinite Bethe roots seem to be absent for all spin-sectors when A = q ±s * , 
i.e. for the boundary conditions where the symmetry generators of the loop algebra 
are known for all spin-sectors [0.. There appears to be a close connection between the 
form of the symmetry algebra generators and the presence of infinite or vanishing Bethe 
roots. This might hold the clue for constructing the symmetry algebra outside the 
commensurate sectors when A = 1. Further investigation is needed to clarify this point. 

6 Conclusions 

In this article we have connected two known alternative methods of solving integrable 
models, the algebraic Bethe ansatz and Baxter's concept of auxiliary matrices. Within 
the community of integrable models the notion of auxiliary matrices has historically 
received less attention than the Bethe ansatz, even though it applies to a wider range 
of models. Unlike the Bethe ansatz auxiliary matrix do not rely on the existence of a 
pseudo-vacuum or spin conservation; see the discussion in Section 3. One possible reason 
that auxiliary matrices have not earlier been investigated in more detail, might be that 
Baxter's construction procedure ^H] leads to a final auxiliary matrix whose algebraic 
structure is not of the simple form ()49|) . Notice also that Baxter's auxiliary matrix for 
A = 1 is limited to the case of even spin-chains, M G 2N, when q is not a root of unity 
and to MmodiV e 2N when q N = 1; see condition (9.8.16) in For a construction 
of Baxter's auxiliary matrices at more general boundary conditions see e.g. |34j . 

In contrast the construction of auxiliary matrices using intertwiners of quantum 
groups yields simpler algebraic expressions, has no restrictions on the length of the 
spin-chain and, most importantly, allows for maintaining the Yang-Baxter relation (|15|). 
It is the last fact which put us in the position to link the formalism of the algebraic 
Bethe ansatz with the Q-operator or auxiliary matrix. Furthermore, the use of repre- 
sentation theory in the construction of Q-operators enables one to make direct contact 
with the underlying algebraic structure of the integrable model, the quantum loop al- 
gebra. The auxiliary matrices discussed in the present article are therefore a key link in 
understanding the relation between the Bethe ansatz and the representation theory of 
quantum groups. 

The importance of this aspect of the Q-matrix is particularly highlighted in the root 
of unity case. As pointed out in [2] the auxiliary matrices constructed at a root of unity 
PP might serve as an efficient tool to analyze the irreducible representations of the loop 
algebra symmetry. In our investigation of the spectrum of the auxiliary matrices via the 
Bethe ansatz as well as the fusion hierarchy we came across the polynomial (cf formula 
dSHD in the text) 




-1) 



M 



,N' 



±1 . 



(103) 



P B (zq- 2k )P B (zq 



-2k-2 
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Recall that the Bethe ansatz equations are sufficient to ensure that the right hand side 
of this identity is a polynomial. This expression has been conjectured to be the classical 
Drinfeld polynomial j2Hj. The latter describes the irreducible representation of the loop 
algebra spanning the degenerate eigenspace of the transfer matrix; see [2j for concrete 
examples. 

What is the benefit of the identification? The symmetry algebra can be used to 
compute the Drinfeld polynomial, i.e. the left hand side of equation IjlOHj) . We can 
then use the above identity to extract the Bethe roots for each degenerate eigenspace 
of the transfer matrix. At the moment this identification still awaits proof as well as its 
extension to the quasi-periodic boundary conditions A = q ±s * investigated in [5]. The 
results of this paper provide a further step towards this direction. 

Away from a root of unity our interest has been to make contact with the findings in 
[Tol IT7j for the Coulomb gas formalism of conformal field theory. This also provided the 
motivation to perform the algebraic Bethe ansatz computation of the spectrum of the Q- 
operators for an arbitrary quantum space in order to accommodate the CFT setup. For 
the XXZ spin-chain we started from the results obtained in to analyze the spectrum 
of the auxiliary matrix and to find the analogue of the functional equations reported 
in P?] for the six-vertex model. The noteworthy difference between the constructions 
in and JH] is the occurrence of free parameters in the representation spanning the 
auxiliary space of the Q-operator; see definition (|23j) . We saw that the introduction of 
such parameters is natural in light of the representation theoretic background (f2T)J) of 
the TQ-equation (H; see also [Hj. 

What is the role of these parameters with hindsight of the spectrum of the six-vertex 
model? Our result for the eigenvalues ()66)) and ([T2*]) showed that the parameter r$ in ([2*5]) 
is needed to break spin- reversal symmetry, while the parameters r^a simply reflect the 
polynomial structure of the eigenvalues of the auxiliary matrix in the spectral variable 
z. One of them is needed to truncate the auxiliary space in order to achieve convergence 
when q is not a root of unity. The freedom in choosing the remaining parameter can 
be used to decompose the eigenvalue ([fill]) into the polynomials (175]) and ([7B]) ; see also 
(J77J). This decomposition has been observed in [T7j in the context of CFT (see equation 
(4.10) therein) and was the starting point for the formulation of a series of functional 
equations all of which we recovered in the case of the XXZ spin-chain; see for example 
equation (|HH|) in the text. As these functional relations have been the starting point for 
connecting the spectrum of the auxiliary matrices with ordinary differential equations 
[55] . it is natural to ask whether this can be also achieved for the XXZ spin-chain. So 
far partial results exist at some roots of unity only |36| . 

Convergence problems did not arise in the root-of-unity case and one might ask 
about the root of unity limit of the auxiliary matrix (|66|). While we did not pursue this 
issue in detail it seems at first sight plausible that in this limit the infinite-dimensional 
representation splits up into an infinite number of finite-dimensional subrepresentations. 
Reducing the trace of the monodromy matrix (|13|) to such a subrepresentation one might 
be lead to the conclusion that one ends up with the auxiliary matrices constructed at a 
root of unity in PQ. This is not true. The full set of auxiliary matrices constructed in pQ 
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does not preserve the total spin S z . Here we only dealt with a subset of the auxiliary 
matrices available at a root of unity. These additional Q-operators also contain free 
parameters but their nature is different from the ones in (|23jl. They reflect the enhanced 
symmetry of the XXZ spin-chain at rational coupling, i.e. q N = ±1. The full scope of 
this symmetry is not accessible via the Bethe ansatz whence further work is required 
to determine the spectrum of all auxiliary matrices in pQ as well as the spectrum inside 
the degenerate eigenspaces; see our discussion in Section 4.1.2. 
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A The spectrum of Q. Proof for ub = 2,3. 

We outline the main steps in proving the conjecture (|4fi|) for the Bethe states (IH7j) with 
ub = 2,3. We will make use of the following well-known commutation relation of the 
Yang-Baxter algebra, 

A X B 2 = ^-B 2 A 1 -^B 1 A 2 , D 1 B 2 =^-B 2 D 1 -^-B 1 D 2 , 

021 &21 012 012 

[Di, A 2 \ = C ^B 2 C X - ^B 1 C 2 , [d, B 2 \ = ^ ( a 2jDi _ AlD2) . 

012 012 012 

Here the indices are a shorthand notation for the dependence of the operators (jOJ and 
the Boltzmann weights (fT2*|) on the respective Bethe roots. For instance, 

bij = b(zi/zj), ... and A x = A(zi), B { = B(zi), ... etc. 

To unburden the formulas we also introduce the symbols 

t = 5i(w/zj) _ P^w/z^^w/zi) t = P k (w/zj) 

kl a k (w/zi) ai + i(w / Zi)a k (w / Zi) ' k a k (w/zi) ' 

The case tib = 2. From the relations pUj) . (J3l"f . p2j) and (j33|) one obtains 
QkkB\B 2 = K\ k K 2 kk BiB 2 Q kk 

H — r^" Qkk+iB 2 Ai + k\ k r 2 k+l BiQ kk+x A 2 — r l k+l -^- Q kk+ iBiA 2 

021 021 

— Qk-ikBiDj. - K\ k rlB 1 Q k _ lk D 2 + Tl^Q k ^ k B x D 2 

0\ 2 012 

+rl +1 rl Q k ^ lk+1 B 2 d + Al k r 2 k+1 r 2 k B x Q k ^ lk+1 C 2 
A 

; 012 

Moving all the 5-operators past the Q-operators except for the first term, this can be 
rewritten as 

1 ^2 

QkkB\B 2 = A kk A kk BiB 2 Q kk + Q kk+x B 2 Ai + -j^- Q kk+X B\A 2 

— Qk-ikB 2 D 1 — ^ Q k ~ikB x D 2 — rl +2 rl +1 jr^- Q kk+2 A X A 2 

- rM-i]^ Qk-2kD x D 2 + ^ti Q A _ lfe+1 A 2 A + %^ Q fc _! fe+ iAA 



+ ^^+iA^\r Qk-ik+iB^ + rj^ +1X r%- Q k _ lk + X B X C 2 

ll kk+l" 12 ^kk+l" 21 

r l+2 r l+i r k^- Qk-ik+2CiA 2 - ri +2 r 2 k+l r 2 k x Akk Q k „ lk + 2 A X C 2 

kk-\-l k — 1 fc+1 



+ rl +1 rlrl_ 17 ^^Q k - 2k+1 CiD 2 + r 2 k+1 r 2 k rl_ lV ^- Qk-2k+\D x C 2 

Ji k-lk il fc-lfc+l 

_ r k+2 r k+l r k r k-l^ Akk Qk-2k+2C\C 2 
lv k-lk+l 
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Acting with this expression on the pseudovacuum all the terms containing C-operators 
vanish. Taking the trace on both sides of the equation we then obtain 



J2Q^BiB 2 \0) H = J2 A lABiB 2 Q kk \0) n 

k k 

+ J2 r l+i Qkk+iB 2 (A 1 /b 21 - D 1 /b 12 ) |0) 



k 

2 



k 



J2 -j^ Qkk + iB l {A 2 /b 12 - D 2 /b 2l ) \0) H 

012 



+ E ^-^Qk-ik+iAxD, + r -^±l Qk _ lk+1 A 2 D 1 ) \0) n 

V & 21 012 



E U^^Qkk^AA. + rlrl.^Q^D.D,) \0) n 

u \ 1V kk+l ^k-lk J 



Employing the Bethe ansatz equations (|55J). 

(01^10)^6^ = 6^ (0|A|0)„, z,j = l,2, 

and using the identity 



12 2 1 



A 1 A 1 



fc+l'fc , ' fc+1' fc _ 1 2 ""-fc— lfc— 1 2 1 - 1 v fc+l _ f> 

i ; r fc+i r fe~rL r fe+i r fc~ri — u 



°12 °21 A fc-lfc A fcfc+1 

we see that the remaining terms cancel. This proves for % = 2. 



The case % = 3. The calculation follows the same strategy as before. The various 
steps can be simplified by observing that the expressions must be symmetric under 
any permutation of the Bethe roots. We omit the details of the computation and only 



30 



provide some intermediate steps. Starting from the result for the case % = 2 one has 



QkkB\B 2 B 3 |0) w = A 1 kk Al k B 1 B 2 QkkB3 \0) n 

+ Qk-ik+i^ {^B 3 A lD2 - ^B^D, - ^t 3 B 2 A lD3 + gfeB^Ds) \0) n 
+ Qk^k + :^T (sfc^W* - ^B 2 MD 1 - ^fe^ 2 £>3 + Sif^A^a) |0> w 

+ Qk - ik+i (4ifc r '+ ir ^ i(A3jD2 - + 4fcfc r ^ ir ^ 2(A3jDi - AijDs) ) |0) * 

- Qkk + 2^^ {^BsAA, - ^B lA3 A 2 + ^B 1 A 2 A a - j^B^A,) \0) n 

- Q^^tr (fcfc^iA, ~ ^t B ^ D * + tt B ^ - iiL B ^ D *) \°)n 

~ G*-i fc+2 ( A g;^^ i(A 3 D 2 - A 2 D 3 ) + If^ ^W - A^s)) |0> w 

( A %^l:f 3 D ^ A ^ - A ^ + A[ t;£ c ^ (^A - ^s)) io>„ 

Here we have dropped all terms where C acts on the pseudovacuum first and those 
which contain more C than 5-operators. The last identity is further simplified to 

QkkB±B 2 B 3 \0) H = A 1 kk Al k Al k B 1 B 2 B 3 Q kk \0) n 

+ Qkk+i (^B 2 B 3Al + p^ Bl B 3 A 2 + ^ Bl B 2 A 3 ) \0) H 
\b 2 ib 3 i o l2 b 32 613623 / 

- Qk-ik (r^-B 2 B 3 D 1 + -^-B ± B 3 D 2 + -^-B 1 B 2 D 3 ) \0) n 
V 613612 6 2 i6 23 631632 / 

-^-S^Ua + -i^SsAaA + ... |0> w 
621631623 b l2 b 32 b 13 J 

~ ^ f ^i fcr H 2 1 +1 ^^i^2 + \0) H ~ Qk-2k { ^ A f} B 3 D lD2 + |0> w 
V A fcfc+i fo 3i632 / \A l k _ lk b l3 b 23 J 

-9^ ( %%^ A 1 ^ 3 + ...) |0> w 

(A 1 />< 1 y< 2 y, 3 \ 

A fc -lfc ft 13023 / 

+ A t f rUA^rt+sQkk+sA^A, \0) n - f kk f k ^.^W |0> w 

7i fcfc+l ii A:fcH-2 IV k-lk^k-2k 

The omitted terms in the parentheses are obtained by symmetrization w.r.t. the Bethe 
roots. Again one computes that upon taking the trace on both sides of the equation 



31 



and invoking the Bethe ansatz equations, 

(oi^io) H n^=(°iAio) w n^' 

all terms vanish except the first one. This proves (J46j) for % = 3. We leave the case 
for general rig to a future calculation. Here we shall be content with supporting the 
conjecture for n# > 3 by making contact with various functional equations for the 
eigenvalues; see (|55 Jl . (foT|). p ty . (ITUUJ) and flTTTTl) . 
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